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Abstract 



K^ We study the thermodynamic stabihty of charged black holes in gauged 

1/-^ ' supergravity theories in Z) = 5, Z? = 4 and D = 7. We find explicitly the 

0^ , location of the Hawking-Page phase transition between charged black holes 

^ , and the pure anti-de Sitter space-time, both in the grand-canonical ensemble, 

^^ I where electric potentials are held fixed, and in the canonical ensemble, where 

total charges are held fixed. We also find the explicit local thermodynamic 
stability constraints for black holes with one non-zero charge. In the grand- 
canonical ensemble, there is in general a region of phase space where neither 
Qh. the anti-de Sitter space-time is dynamically preferred, nor are the charged 

r~| ! black holes thermodynamically stable. But in the canonical ensemble, anti-de 

Sitter space-time is always dynamically preferred in the domain where black 
holes are unstable. 

We demonstrate the equivalence of large R-charged black holes in D = 5, 
5t , D = A and D = 7 with spinning near-extreme D3-, M2- and M5-branes, 

respectively. The mass, the charges and the entropy of such black holes can 
be mapped into the energy above extremality, the angular momenta and the 
entropy of the corresponding branes. We also note a peculiar numerological 
sense in which the grand-canonical stability constraints for large charge black 
holes in D = 4 and D = 7 are dual, and in which the D = 5 constraints are 
self-dual. 



I. INTRODUCTION 

Since the work of [1] on the thermodynamics of near-extreme D3-branes, there has been 
considerable interest in the extent to which supergravity solutions reflect the thermal prop- 
erties of non-abelian gauge theory in 3 + 1 dimensions with A/" = 4 supersymmetry. The 
conformal invariance of the underlying theory suggests that the only part of the supergrav- 
ity solution that is relevant is the near-horizon region, which amounts to a black hole in 
five-dimensional anti-de Sitter space (AdS^) times a five-sphere (S^). Indeed, following the 
conjecture [2] that string theory on AdS^ x S^ is physically identical to A/" = 4 gauge theory 
on the boundary of AdS^, it was shown [3,4] that the thermodynamics of large Schwarzschild 
black holes in AdS^ matches the expected thermodynamics of the gauge theory, at least up 
to constants. In part this amounts to a rephrasing of [1]; but also [3] went on to propose 
a description of pure QCDg in terms of the finite-temperature theory on a near-extreme 
D3-brane. A key ingredient to the confinement-deconfinement transition described in [3] 
was that the gauge theory is on a three-sphere of finite volume. The finite volume breaks 
the conformal invariance which would otherwise make a phase transition impossible. 

In [5], the study of near-extreme D3-brane thermodynamics was extended to include 
the case where the branes have angular momentum in a plane orthogonal to the branes. 
We call this type of angular momentum spin. Like most charged black hole solutions, the 
D3-brane with no spin exhibits positive specific heat sufficiently close to extremality. This 
property has long been regarded as the sine qua non for comparison with field theory, since 
field theories are always expected to have positive susceptibilities. The surprise in [5] was 
that even in the near-horizon limit where field theory is supposed to be applicable, the 
thermodynamics is not always stable. If the spin is too large compared to the energy, the 
near-extreme brane solution no longer represents even a local maximum of the entropy; 
the entropy fails to remain a subadditive function of external variables such a energy and 
angular momentum. The stability analysis has been generalized to the M2- and M5-branes, 
and to multiple angular momenta [6]. For D3-, M2- and M5-branes, Cai and Soh [7] have 
discussed the critical behavior of specific heats and other susceptibilities near the boundary 
of stability, and also at other points in phase space related to stability with fixed angular 
momenta. 

Subadditivity of the entropy is a generalization of the requirement of positive specific 
heat. It amounts to the "effective" Euclidean action for the black hole configuration being a 
local minimum. Thus the stability criterion used in [5] is rather different from [3] , where the 
location of the confinement-deconfinement transition was found by comparing the Euclidean 
action of two topologically distinct geometries which made competing contributions to the 
path integral. One geometry represented a black hole, while the other represented a gas 
of particles in anti-de Sitter space. The discontinuous transition between the two was first 
studied by Hawking and Page [8]. 

One of the goals of the present paper is to compare the Hawking-Page transition to the 
subadditivity condition on the entropy, for charged black holes in anti-de Sitter space, i.e., 
charged black hole configurations in maximally supersymmetric gauged supergravity [9-11]. 
In a realization of the gauged supergravity theories as a truncation of string theory or 
M theory compactified on a sphere, the gauge group is the isometry group of the sphere. 
In the AdS/CFT correspondence this gauge group becomes the R-symmetry group of the 



boundary CFT which is thought to be dual to string theory on the anti-de Sitter geometry. 
(Recall, the R-symmetry algebra is by definition that part of the super symmetry algebra 
which commutes with the Poincare symmetries but not the supersymmetries. It is like a 
bosonic flavor symmetry, but evades the Coleman-Mandula theorem by nevertheless partic- 
ipating in the same superalgebra as the Lorentz group generators). We therefore refer to 
this type of charge as R-charge. Our work extends that of [5] in several ways. First, we 
show that the Kaluza-Klein reduction of near-extreme spinning branes in fact corresponds 
to large R-charged black holes in D = 5, 4 and 7. Their mass and charges correspond to 
the energy above extremality and the angular-momenta of the near-extreme spinning D3-, 
M2- and M5-branes. For general R-charged black holes we find explicitly the constraints 
for the Hawking-Page phase transition as well as subadditivity conditions in the case of a 
single R-charge turned on. We do the analysis with the temperature and the R-charges 
fixed, and also with the temperature and the electric potentials fixed. The electric poten- 
tials are the thermodynamic dual variables to the R-charges. We refer to the first type of 
analysis as stability in the canonical ensemble, and to the second type as stability in the 
grand-canonical ensemble. A subtlety regarding the differences between ensembles and the 
relevance of thermodynamic variables versus fixed parameters is explained in section III A. 

The general conclusion from local stability analyses is that branes which spin too fast, 
or equivalently black holes with excessive R-charges, are unstable. This begs the question, 
what do such configurations break into? If the conserved R-charges are regarded as fixed 
parameters (canonical ensemble), then in view of our results it seems plausible that exces- 
sively R-charged black holes in anti-de Sitter space undergo a Hawking-Page transition to 
a gas of particles. Our analysis shows that this is indeed the case. The transition is an 
essentially quantum phenomenon which eliminates an event horizon and a singularity from 
the space-time. On the other hand, if instead of the total R-charges their dual potentials are 
specified at the boundary of anti-de Sitter space (grand-canonical ensemble), we find that 
the Hawking-Page transition cannot altogether explain the instability of excessively charged 
solutions. We show that for the grand-canonical ensemble there is always a domain of large 
R-charges for which the black hole is dynamically preferred in the sense of Hawking and 
Page, but is locally thermodynamically unstable. That is, the black hole's Euclidean action 
is smaller than that of the gas of particles in anti-de Sitter space, but nevertheless is only a 
saddle point and not a local minimum of the action. 

In Section II we exhibit the R-charged black hole solutions in AdS^, found in [10], which 
will be the focus of the paper. In Section HI we describe the conditions for local thermo- 
dynamic stability and for the Hawking-Page-type transition. In Section IV we demonstrate 
explicitly that spinning near-extreme D3-branes can be Kaluza-Klein reduced to large R- 
charged black holes in AdS^, and show (Subsection IV B) how the thermodynamics of spin- 
ning D3-branes can be recovered from our results in Section III. In Sections V and VI we 
extend our results to charged black holes in AdS^ [11] and AdSi, respectively. In Section VII 
we relate these black holes to spinning M2-branes and M5-branes. 



II. BLACK HOLES IN FIVE DIMENSIONAL GAUGED SUPERGRAVITY 

In this section we summarize the black hole solution in D = 5 N = 8 gauged supergravity. 
This solution was found in [10] as a special case (STU-model) of the solutions oi D = 5 N = 2 
gauged supergravity equations of motion. The relevant bosonic part of the Lagrangian can 
be cast in the following form [10]: 

(1) 

The space-time indices fi, z/, etc. run from to 4. The metric has signature -++++. R is 
the scalar curvature, F* {i = 1, 2, 3) are the Abelian field-strength tensors, e = y/—g is the 
determinant of the Fiinfbein Cm"", and V is the scalar potential given by [10]: 

V = 2E^ (2) 

The metric on the relevant part of the scalar manifold is 
Gij = |diag((X^)~^, (X^)"^, (X^)"^). Here X* are three real scalar fields X\ subject to 
the constraint nf=i-^* = 1- The solution is specified by the following metric: 

ds'' = -{H.H^H.r^/^fdf + {H.H^H^fl^ (j-^dr^ + r^rff^g,^) , (3) 

where 

f = k--^ + g'r^HiH2Hs, H, = 1 + ^, ^ = (1,2,3), (4) 

and (ifig ^ is the metric on S^ with unit radius if /c = 1, or the metric on R^ if A; = 0. 9 = j^ 
is the inverse radius of AdS^, and related to the cosmological constant A = —6g^ = —ji- 
The three real scalar fields X* and three gauge field potentials A^ are of the form 

X^ = Hr\H^H,H,f' , A] = -^ ,^ = 1,2,3 . (5) 

Above, we have chosen the Newton's constant Gs = ^} Here the following parameterization 
for g, (the "charges" entering the metric) and qi (the physical charges) is introduced: 

qi = II sinh^ /3j , qi = fi sinh (3i cosh /3j . (6) 

In the following we shall concentrate on the black hole solutions with k = 1. The solution 
is specified by the ADM mass: 
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where the second equahty follows from the definition of r_|_, the radial coordinate at the 
outer horizon, as the largest non-negative zero of the function /(r) appearing in (4). 

When all the three charges are turned it is not possible to obtain the explicit form of 
the entropy S and the inverse temperature j3 in terms of physical charges qi and the ADM 
mass M . However, it turns out that the entropy S and the inverse Hawking temperature 
(3 can be expressed explicitly in terms of the "charges" qi and the radial coordinate of the 
outer horizon r+: 
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(8) 



i=l 



P=7^ 



2nrlJuUirl + q, 



Th 2r% + ri{l + Y.Uq^)-Y\Uq^ ' 



(9) 



Again, note that qi rather than qi are the physical charges (i.e., the conserved charges to 
which Gauss's law applies). One can express 



q^i =qi{rl + qi 



i+i^n(^++^^-) 



j¥=i 



(10) 



Also, the electric potentials at the outer horizon take the following form: 

Qi 



= Alir^) 



ri + qi 



(11) 



The parameterization (8), (9), (10), (11) will turn out to be useful when deriving the the 
thermodynamics stability constraints and constraints for the Hawking-Page phase transition 
for these configurations. 

In the case of a single non-zero charge, say, qi = q ^ 0, the horizon coordinate r+ and q 
can be expressed explicitly in terms of the ADM mass M and physical charge q, and thus all 
the corresponding thermodynamic quantities S, [3 and can be written explicitly in terms 
of M and q. The following relationships are useful in this case: 



l-g+ J(l + g)2 + 4/i 






g- + ^-^ 



-M 



M2 + 3g2 . 



(12) 



Using these expressions, one can express the entropy S = 27rr^wr^ + q wholly in terms of 
M and q. 



Before we proceed, a comment on the choice of units is in order. It is not completely 
general to choose units where simultaneously G^ = ^ and g = j = 1. (L is related to the 
asymptotic cosmological constant A = —^■) The reason is that G5 is a length cubed and 
g is an inverse length, so G^g^ is a pure number. If A^ is the number of units of Dirac flux 
supporting the AdS^ x S^ geometry, then one finds G^g^ = ^f^. So picking G5 = | and 
g = 1 simultaneously actually fixes A^ = ^/2. 

Nevertheless, in the following sections we shall take G*5 = f and L = - = 1. Powers 
of L can be restored by making the following replacements: r^-^^,qi^j^,fi^j^, 
M — i> ^, and qi ^ jj- In Section IV, where we relate the spinning D3-branes and their 
thermodynamics to that of large black holes, we restore L explicitly. 



III. PHASES OF R-CHARGED BLACK HOLES 

A. Ensembles and the nature of stability 

We will explore two notions of thermodynamic stability: local and global. Roughly 
speaking, the distinction is local versus global minima of an appropriate thermodynamic 
potential. There are various subtleties to explain, and we will address them in this section. 

Local thermodynamic stability is the statement that the entropy S* as a function of the 
other extensive thermodynamic variables Xi is subadditive in a sufficiently small neighbor- 
hood of a given point in the phase space of possible Xj. In our case the variables Xi are M 
and the conserved charges g^. When S' is a smooth function of the Xj, subadditivity is equiv- 
alent to the Hessian matrix [ gf.^ J being negative definite. Again provided S is smooth, 
the region of local stability is bounded by a locus of simple zeroes of the determinant of the 
Hessian (or possibly zeroes of higher odd order — but this is a non-generic situation that we 
will never encounter). 

On can also characterize the region S of phase space in which local thermodynamic 
stability holds in terms of Legendre transforms. S is the largest possible region satisfying 
two criteria. First, it contains a particular reference point which we know is stable: for 
instance, a point where M is very large and all qi = 0. Second, in Legendre transforming 
the thermodynamic potential with respect to some or all of the Xi, the map from Xi to the 
thermodynamic conjugate variables Xi must be one-to-one, so that the Legendre transform 
with respect to the chosen Xi composed with the Legendre transform with respect to the 
conjugate Xi gives the identity (i.e., the Legendre transforms are invertible). The region S 
can be parametrized in terms of Xi or Xj. Thus we see that thermodynamic stability does 
not depend on the choice of ensemble; rather, it is a criterion for being able to Legendre 
transform freely among ensembles. From a mathematical point of view, different ensembles 
are different methods of calculating the same thing. 

There is one big caveat to the foregoing remarks: we have to decide which quantities 
are thermodynamic variables which could in principle be allowed to vary in an experiment, 
and which are fixed parameters of the system which cannot vary. Striking out entries in the 
list of thermodynamic variables Xi actually decreases the dimension of the phase space. For 
example, if we regarded all the charges g, as fixed parameters, then the phase space would 
be parametrized just by the mass M. However, in practice it is convenient to speak loosely 
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of "phase space" as the space of possible values for all the Xi, regardless of which of them 
we regard as thermodynamic variables. The more Xi we regard as fixed parameters, the 
larger the region of stability. The reason is that we are decreasing the family of Legendre 
transforms that have to be invertible. Fewer conditions means a larger region on which they 
are all satisfied. 

As a form of rhetorical shorthand, we shall indicate which of the Xi we regard as thermo- 
dynamic variables by specifying a particular ensemble. For instance, the canonical ensemble 
is the one obtained from the microcanonical ensemble by trading entropy, or equivalently 
mass, for temperature via a Legendre transform. Thus any reference to the canonical ensem- 
ble really means that we are thinking of the mass as the only thermodynamic variable in the 
microcanonical ensemble, and that the charges g, are fixed parameters. The grand-canonical 
ensemble is obtained from the microcanonical ensemble by Legendre transforming with re- 
spect to entropy and all the charges, so we are now thinking of the charges as thermodynamic 
variables, too. 

The local stability computation in any ensemble — where by speaking of one ensemble or 
another we are distinguishing thermodynamic variables from fixed parameters, as explained 
in the previous paragraph — can in principle be carried out by finding zeroes of the determi- 
nant of a submatrix of [ gf,£.. ]- The appropriate submatrix is the one which includes only 
those rows and columns pertaining to the Xi which are thermodynamic variables. For the 
canonical ensemble, this determinant is essentially the specific heat at constant qi. For the 
grand-canonical ensemble, it is a combination of this specific heat with other susceptibili- 
ties. Direct computation of these quantities tend to be extremely burdensome. We find it 
more efficient to work with the thermodynamic potential of the specified ensemble, but as a 
function of the extensive thermodynamic variables rather than the thermodynamic variables 
which are usually regarded as the independent variables for that ensemble. The method is 
best illustrated by an example. Suppose we are working in the canonical ensemble, so the q^ 
will always be regarded as fixed (and from here on notationally suppressed). Consider the 
Helmholtz (Euclidean) action: 

Ih = I3F = I3M-S . (13) 

Suppose we write M explicitly in terms of S: this amounts to a complete specification of 
the thermodynamics in the microcanonical ensemble. Thus we write Ih as a function of S 
and the free parameter /3, which will be determined by the Legendre transform condition: 

If we take a second derivative, 

[dS^ J^ ' dS^ {dS/dMf ' ^ '' 

(where we have used (14) in the second equality), then we see that (^|f-) = is equivalent 
to ^^p- = 0, which is the relevant local thermodynamic stability equation. In fact, (15) is 
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still rather tedious to compute, but a change of variables makes life much easier: we find we 
can write S wholly in terms of r_|_, which is the location of the outer horizon (see equation 
(8)). It does not matter that the g^ enter into the relation between S and r+, since they are 
just fixed parameters. Next we note that the condition 

in the presence of the constraint (14) is equivalent to the vanishing of (15), provided dS/dr^ 
is non-singular. Equation (16) represents the essence of the method we use to evaluate local 
thermodynamic stability. Working in the grand-canonical ensemble, there are more thermo- 
dynamic variables {qi as well as the mass or the entropy), and the convenient reparametriza- 
tion is in terms of r+ and qi. Working through the steps (14), (15), and (16), it is straightfor- 
ward to show that the zeroes of the determinant of the Hessian of S with respect to M and 
the (ji coincide with the zeroes of the determinant of the Hessian of the Gibbs (Euclidean) 
action Iq, written in the form 

lG = P{M-J2q,<Pi)-S , (17) 

with respect to r^ and g^. Where -^ appears in (16), it will be replaced in the grand- 
canonical analysis by the Jacobian J' = det g/ ''^^K . It is straightforward to check that this 
Jacobian is non-singular. This completes our discussion of local thermodynamic stability. 

The Euclidean version of an R-charged black hole in AdS^ is characterized by some 
boundary data: the period for Euclidean time, and the potential, which in ten-dimensional 
language amounts to the angular velocity with which the S^ is rotating. There is another ge- 
ometry with the same boundary data: perfect, Weyl-flat, Euclidean AdS^ times an S^ which 
again is rotating with the angular velocity specified by the potential. In five- dimensional 
terms, this geometry is interpreted as a gas of R-charged particles in anti-de Sitter space. 
Following Hawking and Page [8], we can ask whether this geometry or the black hole geome- 
try is favored. The geometry with smaller action will dominate the Euclidean path integral. 
The nature of the boundary data determine how the action must be computed. If it is indeed 
the potential and the temperature which are specified at the boundary, then the relevant 
action is the Gibbs action, (18). If instead the total R-charge is fixed, it is appropriate 
to compute the Helmholtz action, (13). We will discuss the computation of these actions 
in more detail in the next two Subsections. In both cases our strategy will be to subtract 
off the pure AdS^ action from the black hole action to obtain a finite quantity, Iq or In- 
The vanishing of these actions will indicate the Hawking-Page phase transition. Past this 
transition point the black hole cannot be the global minimum of the Euclidean action. If 
the pure anti-de Sitter geometry is the only geometry which competes with the black hole in 
the path integral, and if the black hole is locally stable, then the Hawking-Page transition 
represents the edge of the region where the black hole is globally stable. 

To summarize, in judging local stability the choice of ensemble amounts to a decision 
regarding which quantities in the system are thermodynamic variables which are allowed to 
vary within the system, and with respect to which Legendre transforms should be invertible; 



whereas in locating the Hawking- Page transition, the choice of ensemble is motivated by the 
more standard notion of whether a given thermodynamic quantity or its conjugate (in the 
sense of Legendre transforms) is controlled from outside. 

It is a matter of physical motivation which ensemble one decides to use.^ For instance, 
in the study of spinning branes [6] where the world-volume is infinite, any given part of 
the brane can be regarded as a "system" in the presence of a "heat bath" which is the rest 
of the brane. The boundary between the system and the heat bath is fictitious, and both 
energy and charge can flow across it. This suggests that for the purpose of judging local 
stability, one had better use the grand canonical ensemble. Using the canonical ensemble, 
where Legendre transforms with respect to the charge density are disallowed, one is ignoring 
the possibility that the charge density might become inhomogeneous. 

Holography also seems to favor the the grand-canonical ensemble, since the temperature 
and electric potentials (which are fixed in the grand-canonical ensemble) can be specified in 
terms of the asymptotic geometry. Specifically, the inverse temperature is the circumference 
of the Euclidean time's S^, and the electric potentials divided by the temperature are the 
twist one makes before identifying around this circle. The standard AdS/CFT approach 
would be to sum over all geometries with the given asymptotic circumference and twist. 
Assuming that only the black hole geometry and the pure anti-de Sitter geometry are com- 
peting in the sum, one judges which is preferred according to which has a lower Gibbs free 
energy. From this point of view the grand-canonical ensemble is the relevant one, and we 
study it in Subsection HI B by computing the Gibbs Euclidean action. Partly to illustrate 
some interesting differences between the ensembles, we calculate the Helmholtz Euclidean 
action in Subsection IIIC and discuss the stability constraints that follow from it. 

B. Grand-Canonical Ensemble 

In this subsection we spell out the properties of the Euclidean effective action for the 
grand-canonical ensemble of R-charged black holes discussed in Section II. A state in the 
grand-canonical ensemble is usually specified by the inverse temperature f3 and the electric 
potentials Al = (pi {i = 1,2,3). 

The Gibbs Euclidean action can be obtained following the procedure, initiated by Hawk- 
ing and Page [8] for the Schwarzschild black hole in ^^5*4, further developed by York et 
al. [12], and applied to D = 4 anti-de Sitter Reissner Nordstrom black holes in [13] (see also 
[14]). When employing this procedure the boundary conditions (including the local temper- 
ature) are defined at a finite value of the radial coordinate r^; these boundary conditions 
in turn uniquely determine the reduced Euclidean action which has both bulk and surface 
term contributions. One also subtracts from it the pure AdS^ space-time contribution. As 
the last step the limit r^ ^ oo is taken, i.e., the artificial boundary is removed. For further 



^In the context of spinning branes the distinction between the canonical and grand-canonical 
ensembles and the correspondingly different local stability constraints were addressed by Cai and 
Soh [7]. 



details see, e.g., [12,14,13].^ 

The above procedure will assign zero action to the pure AdS^ geometry. With no sub- 
tractions and no cutoff, the action integral formally diverges, so some choice of zero point is 
necessary. The AdS^ geometry also involves a constant electric potential. Constant electric 
potential yields zero field strength and hence makes no contribution to the action, and the 
the anti-de Sitter geometry is undistorted. In the ten-dimensional geometry, the constant 
electric potential is interpreted as a twist on the S^ part of the geometry. To understand 
this, note that the gauge fields arise from mixed components of the ten-dimensional metric, 
Qfj^a, where n is an AdS^ index and a is an S^ index. For instance, if we turn on a single con- 
stant electric potential, corresponding to rotations in the direction of an S^ parametrized 
so that the metric is rfs^ = dO"^ + siv? Odcj)^ + cos^ 6 dVL\, then what we are doing from a 
ten-dimensional perspective is making gtff, non-zero. This off-diagonal entry in the metric 
can be exactly killed if one changes the angular variable from to = — Vlt. Vl has the 
interpretation of an angular velocity at which the S^ rotates. Now, is single valued in the 
periodized Euclidean geometry, whereas jumps by ^/Th at the identification point. This 
is what we mean by twist on the S^. If one applies Kaluza-Klein reduction to get back to 
five-dimensions, then in the normalization conventions used in Section IV B, we have the 
electric potential Af = — i7. 

The Gibbs Euclidean action takes the following form: 

IG = (3{M-J2Q^<P^)-S (18) 

where M is the ADM mass, qi is the physical charge and S is the entropy. S = j^ where 
A is the area at the outer horizon. 

The structure of this action is expected from thermodynamic considerations as well. The 
physical, Lorentzian, black hole solutions are extrema of (18). The expressions (9) and (11) 
for the inverse Hawking temperature (3 and the electric potentials 0i at the outer horizon 
can be obtained by extremizing (18) with respect to S and qi with 0j and (3 held fixed. As 
discussed at the beginning of this Section, these extrema can also be obtained by finding 
zeroes of the first derivatives of (18) with respect to qi and r^, again keeping (3 and 0j fixed, 
provided the Jacobian J' = det g/ ''^'l is non-singular. 

For the R-charged black holes of Section II the Gibbs action (18) can be expressed 
explicitly in terms of q^ and r_|_, and it takes the following suggestive form: 



^In Ref. [3,15] an alternative approach in the study of the of AdS-Schwarzschild black holes, 
suitable for the study of the higher loop (a') corrections to the effective action, was developed. 
There one considers only a bulk term, but again there is a cutoff at large radius both for the pure 
AdS geometry side and for the AdS-Schwarzschild black hole side. The pure AdS geometry is 
periodized with a temperature in such a way that at the cut-off radius the circumference of the 
compactified S^ is the same as for AdS-Schwarzschild black hole. Then by subtracting the pure 
AdS bulk integral from the AdS-Schwarzschild black hole bulk integral, one obtains a finite result 
corresponding to the AdS-Schwarzschild Euclidean action. 
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^^' = Ak-mri+J , (19) 

+ L i=l 

where P is determined by (9). (We have set L = 1, but powers of L can be restored using 
the replacements described at the end of section II). 

Now let us specialize to only one non-zero charge: say Q'l = Q' 7^ but q2 = Qs = 0. 
The local thermodynamic stability is most efficiently analyzed by computing the Hessian of 
(18) with respect to r+ and q with f3 and held fixed, in analogy with (16). The stability 
constraint is 

2rl + rl{q + l)-{q-lf>0. (20) 

In the phase space parametrized by (M, g), or more conveniently by {r^,q), the zeroes of 
the left hand side of (20) determine the critical lines which form the boundary of stability. 
There are two branches: 



g±(r+) = l + lrl± iy'Sr^ + Orj . (21) 

The regions for which q-{r^) < q < q'+(r_|_) correspond to locally stable regions, i.e., Iq^ is 
a local minimum. When g_ = the critical point is at r+ = 4^ [q = 0, M = |] [10]. The 
critical lines merge at g+ = g_ = 1 and r_|_ = [M = q = 1]. This point is a BPS solution 
with zero entropy and finite (3 = n. Black holes in the regions q < g_(r+) and q > g_|_(r_|_) 
are thermodynamically unstable. The region q < g„(r_|_) charged black holes which are 
much smaller than the size of anti-de Sitter space: r_|_ <^ 1, g <^ 1. The curvature of AdS 
is negligible in the vicinity of these black holes, so we are predicting that singly charged 
black holes should be unstable even in fiat space. And they are: they can be obtained in a 
T^ compactification of type IIB string theory as NS5-branes or D5-branes wrapped around 
the T^. The thermodynamic instability of NS5-branes is familiar: the temperature in the 
limit of extremality is the Hagedorn temperature of the fractionated instanton strings on 
the world- volume [16], and as non-extremality increases the temperature goes down. 

The location of the Hawking-Page phase transition between the black hole solution and 
the pure AdS^ solution is determined by the constraint Iq^ = 0: 

go(r+) = l-r^, (22) 

where we have used (19). For q > q'o(^+) the black hole solution has Iq^ < and thus it is 
the dynamically preferred solution. go(^+) = corresponds to r+ = 1 [q = 0, M = 3]. This 
is the Hawking-Page phase transition for the D = 5 ArfS'-Schwarzschild black hole solution. 

For q > q'+(r+) > qo{r+), black holes are favored over the pure anti-de Sitter space-time, 
but they are locally unstable. In this case our evaluation is that we have failed to find 
the entropically preferred state of the system. One possibility that suggests itself is that a 
multi-black hole solution might be favored, but it seems difficult to imagine a static solution 
describing multiple black holes in AdS^. 

For g_(r+) < q < qQ{r+), pure anti-de Sitter space is entropically favored over the black 
hole solution, but the black hole solution still corresponds to a local minimum of the Gibbs 
action. In the uncharged case, this is the region between r+ = -^ [M = |] and r+ = 1 
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FIG. 1. In Figure la stability domains are plotted as q vs. r_(_ for the grand-canonical ensemble 
of D = 5 R-charged black holes with one charge turned on. In Figure lb the same stability plots are 
exhibited, but now as q (the physical charge) vs. M (the ADM mass). Only the region satisfying 
the BPS bound, M > q, corresponds to physical black hole solutions. The vertically shaded areas 
correspond to the regions where the AdS^ is the preferred solution and the horizontally shaded area 
correspond to the regions where the black hole solutions are local minimum of the Gibbs action. 
The checkered area is the domain of common overlap: there black holes still correspond to the 
local maxima of the entropy (local minima of the Euclidean action); however, anti-de Sitter space 
is preferred globally. In the unshaded area, neither the black hole solution nor pure anti-de Sitter 
space are global minima of the action. 

[M = 3]. These black holes could be termed "meta-stable" in the sense that they can exist 
in equilibrium with a heat bath of particles in anti-de Sitter space and are stable under small 
perturbations, but the system as a whole could gain entropy by converting the entire black 
hole into a gas of R-charged particles. 

In Figure la we show all these stability domains in a plot of q vs. r+. To help the 
reader's intuition we also exhibit the same stability plots in Figure lb, but now as q vs. M. 
(On this graph only the region M > g is consistent with the BPS bound.) The vertically 
shaded areas correspond to the regions where the AdS^ is the preferred solution, and the 
horizontally shaded area correspond to the regions where the black hole solutions are local 
minimum of the Gibbs action. 

C. Canonical Ensemble 

We now repeat the stability analysis for the canonical ensemble; along with the fixed 
inverse temperature [3 , one now fixes the physical charges g^, i.e., the field strengths F^^, 
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(and not the electric potentials (pi) are fixed. Following a procedure analogous to the one 
spelled out for the grand-canonical ensemble in Subsection IIIB, but now with the new 
boundary conditions, one obtains the following form of the Helmholtz Euclidean action: 

lH = f3M-S , (23) 

where M and S are the mass and the entropy of the system and (3 is the inverse temperature. 
Extremizing the above action with respect to S, or equivalently, with respect to r+ (while 
keeping the physical charges qi fixed) determines the correct inverse temperature (3 of the 
black hole solution (9). 

To locate the Hawking-Page transition, we write the Helmholtz action (23) for the black 
hole solution described in Section II as a function of qi and r+: 



/ 



BH P 



^ "2H 



3 3 3 

4/1 „2 \ , „2 



r+ 



l_rl)+rlY,q,{2 + rl) + ^rlY,q,qj+b\{q, 

i=l i<j i=l 



(24) 



where (3 is defined in (9). The phase transition between the black hole solution and the pure 
AdS^ solution is determined by the zero of I^^ (24): 

^o(r.) = '-^f^ . (25) 

In the domain q > qo{r+) the black hole solution has Ih < and thus it is the dynamically 
preferred solution. (For qo = the critical point is again r_|_ = 1 [q = 0, M = 3].) 

To evaluate local thermodynamic stability, we use the method outlined in Subsec- 
tion III A around equation (16). In the case of a single non-zero charge, say qi = q y^ 0, the 
local stability constraint leads to: 

2rl + rl{5q+l)-q'^ + 6q-l>0 . (26) 



This case is simple enough that one can easily verify (26) by explicitly computing [-§j^ 
The zeroes of the left hand side of (26) form two critical lines: 



Q 



g± (r+) = 1 + Irl ±^32 + 64^ + 33r| , (27) 

and the stable region corresponds to the domain g_(r+) < q < q+{r+). Again, g_ = 
corresponds to the critical point: r+ < 4= [g = 0, M = |] [10]. On the other hand now the 

stable domain with r_|_ = corresponds to the range of charges: 3 — 2\^ < q < 3 + 2\f2. 
[Equivalently, this is the region with M = q (BPS-limit) and 3 — 2v^ < g = g < 3 + 2\f2\ 
Black holes in the region q < g_(r_|_) and q > g+(r+) are not thermodynamically stable. 

In Figure 2a we show the stability regions in a plot of q vs. r_|_. In Figure 2b we exhibit 
the same stability regions plotting q vs. M. The vertically shaded areas correspond to the 
regions where the AdSr, is the preferred solution and the horizontally to the regions where 
the black hole solution is the local minimum of the Helmholtz action. 

Note three salient differences between the results for the canonical and grand-canonical 
ensembles. First, as expected from the general considerations of Subsection HI A, the region 

13 



q 4-: 





r^ M 

Figure 2a Figure 2b 

FIG. 2. In Figure 2a we plot the stability domains as q vs. r_|_ for the canonical ensemble of 
D = 5 R-charged black holes with one non-zero charge. In Figure 2b the same stability plots are 
exhibited, but now plotting q vs. M. The vertically shaded areas are the regions where anti-de 
Sitter space is the preferred solution. The horizontally shaded areas are the regions where the 
black hole solution is a local minimum of the Helmholtz action. The checkered area is the domain 
of common overlap: meta-stable black holes in the intepretation of section IIIB. 

of local stability for the canonical ensemble is larger than for the grand-canonical ensemble. 
Second, the region favored by the pure AdS^ is significantly increased for the canonical 
ensemble, and it encompasses the whole region with r+ < 1 and non-zero g^'s. Third, and 
most interesting, the pure AdS^ space-time is dynamically favored in the whole domain of 
the black hole parameters for which the black holes are locally unstable. 

IV. SPINNING D3-BRANES AS R-CHARGED BLACK HOLES 

In this Section we establish the precise connection between the near-extreme spinning 
D3-branes solution and the solutions oi N = 8, D = 5 gauged supergravity found in [10] 
and discussed in Section II. In this Section we restore all factors of L explicitly, but still set 

n 

5—4- 

A. Equivalence of A; = R-charged black holes and near-extreme spinning D3-branes 

Angular momenta in planes perpendicular to the world-volume of D3-branes (which we 
refer to as spins to distinguish them from angular momenta in other planes) are realized 
in the world- volume theory as charges under the global 5*0 (6) R-symmetry group. In the 
reduction of type IIB supergravity on S^, the 5*0(6) isometry becomes the non-abelian gauge 
symmetry of the resulting A/" = 8 gauged supergravity in five dimensions, and the spins 
become 5*0 (6) gauge charges. In the AdS/CFT correspondence, the 5*0 (6) supergravity 
gauge fields in the bulk of AdS^ couple to the gauge theory R-currents on the boundary. 
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It is well-known that the near-horizon geometry of near-extreme D3-branes without spin 
reduces to a Schwarzschild black hole in AdS^. The purpose of this section is to show how 
spin on the D3-branes reduces to charge on the black holes. By explicitly performing the 
Kaluza-Klein reduction on S^, we will obtain the precise relationship between the parameters 
of the R-charged black holes and those of the spinning D3-branes. For the sake of simplicity 
only one angular momentum, so that the black hole will have Q'l = 5' 7^ but q2 = Qs = 0. 
One can verify that the computation extends to the case of all three angular momenta 
non-zero. 

The near-horizon form of the ten-dimensional D3-brane solution with a single angular 
momentum is [17,18] 



Mj^N 



ds^Q = cjMNdx dx 



1 
77 



2£LV, 



2^2 



- hdt^ - """^ '° sin^ edtd(f) + /Ar^ sin^ 
r**A 

+ tdr"^ + dx^ + fAr^dO^ + fr^ cos^ Od^l 
h 



f 



h 



r^A 



A = l + 



' cos^ 9 






A = l + - 



r^A 



h^l(i,'l-< 



A 



(28) 






^^ 



The basic Kaluza-Klein ansatz is [19] 

9mn = 



e'^^^'g^u + A;<A^ 

A 




(29) 



where indices a, /3, 7 run over the compact internal manifold K^ and /x, v run over the non- 
compact spacetime M. The metric on K can vary with the position on M, and the equation 



vol 



K 



j^^vol 



X,oo 



(30) 



defines the "dilaton" field D (not to be confused with the ten-dimensional dilaton 0) in 
terms of how the volume form on K compares to its limit at spatial infinity. The power of 
e^ in (29) is chosen to put the metric g^y in Einstein frame. For the solution (28) we have 



ds^^ = Gnffdx'^dx^ 



„2D 



A 



Ar^de^ + Ar^ sin^ 



' + r^cos^9dnl 



(31) 



A3/4 ■ 

The only non-zero component of A^''' is 



er^ 



LV2A 



(32) 



Finally, the five-dimensional metric g^i, is just the solution (3) with k = 0, fi = r^/ L"^, 
qi = (?, and q'2 = ?3 = 0. Note that although D and ds]^ vary over both M and K, the 
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five-dimensional metric ds\^ and the five-dimensional gauge potential At'^ are independent 
of the internal coordinates. The dependence of D on ii' indicates that the compactification 
is a warped product. 

There is a slight subtlety in comparing the gauge fields in (42) of [10] with (32): super- 
ficially in [10] the field strength seems to vanish because of an explicit prefactor of vk. At 
this point it is useful to recall the approach of [3] where the A; = solution is recovered by 
focusing in on a small, nearly fiat part of the S^ in the k = 1 solution. This amounts to 
taking the limit k -^ 0+. In this limit, q = (? and /x = r^/ L"^ are held fixed. From equation 
(43) of [10] one obtains sinh^/5 = ki'^L'^/r-Q. The factor of k in this relation conspires to 
cancel the one explicitly shown in (42) of [10], and the end result for the gauge field A^ 
pertaining to the charge gi is 

A] = ^A^ . (33) 

The factor of ^ comes from a different choice of normalization of A^ in [10]. For us the 
gauge kinetic term picks up a factor of L from \/GG~^G^^. 

There is some debate over the status of five-dimensional A/" = 8 gauged supergravity 
as a truncation of the Kaluza-Klein reduction of type IIB supergravity/string theory. In 
particular, is it true that solutions of the former can always be promoted to solutions of 
the latter? This seems guaranteed for supersymmetric solutions from closure properties of 
the supersymmetry algebra, but for non-supersymmetric solutions one is making a non- 
trivial statement regarding the equations of motion.'^ Modulo the small technical issue of 
compactifying R^ to S^, non-extreme R-charged black holes provide an example of a non- 
supersymmetric solution of gauged supergravity which can indeed be promoted to a solution 
in ten dimensions: namely the spinning D3-brane solution. The key point is that the gauge 
potential At"^ and the metric g^i, resulting from the Kaluza-Klein compactification of the 
spinning D3-brane solution are independent of the S^ coordinates. 



B. Thermodynamic equivalence of large R-charged black holes and spinning 

D3-branes 

Another way to view the k ^ 0~^ limit is that black hole solutions with A; = can be 
recovered from large black hole solutions with /c = 1 by restricting attention to a angular 
region of the S^ in the k = 1 region where the curvature can be neglected. This was 
explained precisely for Schwarzschild black holes in [3], and the same story goes through 
with charged black holes. New physics shows up in the k = 1 solutions (for instance, the 
confinement-deconfinement transition) because of the finite-size effects of the S^. When 
we take /c — *> O"*", we are making the volume of the S^ very large compared to the cube 
of the inverse temperature. Neglecting all finite size effects should lead us back to the 
thermodynamics of spinning D3-branes whose world- volume is fiat and infinite. The purpose 



^We thank J. Distler for explaining this to us. 
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of this Subsection is to show this exphcitly and to estabhsh the precise connection between 
the respective physical parameters. 

Large R-charged black holes are the ones with r+ ^ L and qi^ L^ . In this limit, the 
expressions for the ADM mass M (7) and the physical charges qi (10) reduce to the following 
form: 

M~|/i, g^~g./x, (^ = 1,2,3), (34) 

where /i is related to r_|_ and g^'s by f{r+) = (see eq.(4)). In the large charge limit this 
equation implies: 

^l[i'■l + q^)-^^irlr^O. (35) 

1=1 

To obtain (35), one can start with /(r_|_)r^ = and drop the term kr^. This makes sense 
in view of the equivalence of /c ^ 0+ and large k = 1 black holes. 
By introducing new variables 

yH = a'^^/llr+, t/i = a^^g^ , (i = 1, 2, 3) , (36) 

where a = (/i^L^)^'^, one can rewrite (35) in the following form: 

n(y| + !/•)- 1/1 = 0. (37) 

This is precisely the "horizon" equation for spinning D3-branes displayed in Section II of 
Ref. [6]. 

One can introduce the energy density e, entropy density s and the charge density ji in 
the following way: 

M 
e = 



Ji = 



Qi _ _£_yi_ 



s=j^ = l^ie^'I'vu , (38) 



where 



\/ = 27^2L^ iV2 = 2L3, 7 = ^^^ • (39) 

The thermodynamic variables introduced in (38), together with the definition of the volume 
V and the flux A^ (39) in terms of L (and G^ if we hadn't set it equal to 7r/4) and also the 
relationship between yu and yi (36), provides a precise translation of the thermodynamic 
variables of large black holes into those of spinning D3-branes. It may seem strange to assign 
a definite volume to the S'^ on the boundary of AdS^. However, as commented on in [20] and 
[21], a choice of radial coordinate is equivalent to choosing a definite metric among those 
with a specified conformal structure. Our choice of radial coordinate was fixed in (3). 
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Even in view of the above translation of variables, the equivalence between the local ther- 
modynamic stability constraints for spinning branes and R-charged black holes is slightly 
non-trivial: equations (34) and (38) define a non-linear relation between e and ji for the 
spinning brane and M and g^ for the R-charged black hole. To see that the stability con- 
straints must turn out the same, think of this relation as a reparametrization of phase 
space, similar to the convenient (r+,gj) parametrization we used in earlier sections. The 
argument in Section III A should guarantee that the stability region does not depend on the 
reparametrization. In the case of one non-zero charge, the stability constraints for R-charged 
black holes for the grand-canonical ensemble (20) and canonical ensemble (26) reduce to the 
following respective forms in the r+/L -^ oo, q/L'^ -^ oo limit: 

2ri — q = -^yff{2 — x'^) > grand-canonical (40) 

2^4 _^ g^2^^ - q^ = ■^yjj{2 + 5a;^ ~ x^) > canonical. (41) 

Here x = ^-. These stability constraints are in precise agreement [5,7] with those of near- 
extreme spinning D3-branes with only one angular momentum turned on. 



V. PHASES OF R-CHARGED BLACK HOLES IN D = A 

So far we have focused our attention on black holes in AdS^ because they are dual to 
thermal states of A/" = 4 supersymmetric gauge theory in 3+1 dimensions. Our analysis is 
extended in this section and the next to black holes in maximally supersymmetric gauged 
supergravity in D = 4 and D = 7. In Section VII we outline their relationship with spinning 
M2- and M5-branes. 

Black holes in D = 4 A^ = 8 gauged supergravity have been studied in Ref. [11]. Their 
Einstein frame metric is of the following form [11,22]: 



ds'^ 



;n H^V'fde + (n H^f/^ (f-'dr' + r'dn,,,) , (42) 



where 



^^ , 2 2 TT TT TT 1 , ^* 



f = k-^+g'r'llH,, H, = l+^-^,{t = l,...4). (43) 



r "^ r 



Here we have chosen 6*4 = |. (Again we shall concentrate on black hole solutions with 
k = 1.) With this choice of G4, g^'s and the physical charges g^'s are defined in (6) and the 
ADM mass is of the following form: 



4 



M = 2/i + ^g, = 2(r+ + rlg^Y[H,) + ^g, , (44) 



i=l j=l i=l 



where in the second equality the radial coordinate at the outer horizon r+ is defined as the 
largest non-negative zero of /(r) (43). 



Again the parameterization of the thermodynamic quantities in terms of qi and r_|_ (in- 
stead of qi and M) is most suitable: 



S 



AGa 



A-n 



\ 






(45) 



/9 



1 



4vrr+Jni=i(r+ + gi 



Th 3rt + 24 Eti g. + n(Ef<, m, + 1) - nti ?. 

The physical charges are then determined in terms of qi and r_|_ in the following way: 



(46) 



;;2 



Qi = <li[r+ + qi) 



1 



1+— n(^++?j) 



jV* 



(47) 



and the electric potentials at the outer horizon take the following form: 



Alir^) 



r+ + qi 



.4). 



(48) 



In the following, along with G4 = ^ we shall also take L 



1. (Note L is related to 



the asymptotic cosmological constant A = — -^0 However, L can be restored by replacing 
r_|_ by ^ and g^ by ^. 



A. Grand- Canonical Ensemble and Stability Constraints 



The Gibbs Euclidean action is of the form: 



IG = P{M-J2Q^<t>^)-S 



(49) 



1=1 



where M is the mass, qi is the physical charge and S = -^^ is the entropy. The extrema of 
this action correspond to the black hole solutions for which the inverse Hawking temperature 
and the electric potentials at the outer horizon are given in (46) and (48), respectively. 

The critical hypersurfaces of thermodynamic stability can again be determined by evalu- 
ating the zeroes of the determinant of the Hessian of second derivatives of (49) with respect 
to r+ and qi. 

The phase transition between the black hole solutions and pure AdS4^ occurs when (49) 
for the black hole solution ceases to be negative. When expressed explicitly in terms of q^ 
and r+, (49) takes the following suggestive form: 



jBH 



/3 



\{{r+ + qi) 



i=l 



(50) 



where j3 is determined by (46). 
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For one non-zero charge only, say gi = g 7^ 0, the the region of local thermodynamic 
stability is specified by 

3r^ + qrl - 2gV^ + 2r^ + 3gr+ - 1 > . (51) 

The zeroes of the above expression correspond to two critical lines: 



f + T^lM 

4r, 4 4 V rj 



g±(r+) = — + ^ ± -J- + 22 + 25r^ , (52) 



with the stable domain satisfying q'-(r+) < q < g_|_(r+). g_ = has a critical point at 
r+ = 4= [8]. Note also that as r_|_ ^ the local stability region is pushed to large q: 
g-t ^ 00 as r_|_ ^ 0. 

The phase transition between the black hole solution and the pure AdS4^ solution is now 
determined by: 

qo{r+) = r+ . (53) 

r+ 

For q > go(^+) the black holes are dynamically preferred over AdS4. As r_|_ -^ 0, the critical 
line is pushed to qo -^ 00. qo = has the Hawking-Page transition at r+ = 1 [8]. 

For q > g+(r_|_) > go(^+) black holes unstable, but at the same time the pure AdS^ space- 
time is not dynamically preferred either. The region g_(r+) < q < go(^+) corresponds to 
the domain where AdS^ is entropically favored, but the black hole solution still corresponds 
to a local minimum of the Gibbs action (49). 

In Figure 3 we plot all these stability regions for q vs. r+. The vertically shaded area 
corresponds to the regions where the AdS4 is the preferred solution and the horizontally 
shaded area corresponds to the region where the black hole solution is the local minimum 
of the Gibbs action. 



B. Canonical Ensemble and Stability Constraints 

Repeating the same analysis now for the canonical ensemble yields the Helmholtz action 
(23) which for D = 4 black holes takes the following form: 



Ih 



BH P 



r I 



i<j i=l i<j<k 4=1 



(54) 



where (3 is defined in (46). I^^ (54) differs significantly from Iq'^ (50). 

In the case of a single non-zero charge the local stability constraint leads to: 

3rl + 6qrl + 2rl + 7gr^ + 2gV+ - r+ - g > . (55) 

The critical line now has only one branch: 
1 



g_(r+) = — i-QrX - ?< + 1 + ,j3Qrl + 60r^ + 21r| - 6r^ + 1 ) , (56) 
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Figure 3 



FIG. 3. The stability domains as q vs. r_|_ for the grand-canonical ensemble of D = 4 R-charged 
black holes with one non-zero charge are given. The vertically shaded areas correspond to the 
regions where the yldS'4 is the preferred solution and the horizontally shaded area correspond to 
the regions where the black hole solutions are local minimum of the Gibbs action. The checkered 
area is the domain of common overlap. The unshaded area corresponds to black holes that are 
unstable but nevertheless favored over pure anti-de Sitter space. Explicit powers of L can be 
restored by replacing r+ with ^ and q with ^. 




Figure 4 
FIG. 4. The stability domains are plotted for q vs. r_|_ for the canonical ensemble of D = 4 

R-charged black holes with one non-zero charge. The vertically shaded areas indicates the regions 

where AdS4, is the preferred solution, and the horizontally shaded area shows the regions where the 

black hole solutions are local minima of the Helmholtz action. The checkered area is the domain 

of common overlap. Explicit powers of L can be restored by replacing r+ with ^, and q with j^. 
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with the region q > q^{rj^) stable. 

The phase transition between the black hole solution and the pure AdS^ solution now 
takes place along the line 

qo{r+) = r\-r+ . (57) 

For q > qo{r^) the black hole solution is the dynamically preferred solution. Pure AdS^ 
is dynamically favored in the whole domain where the black hole is thermo dynamically 
unstable. 

In Figure 4 we plot the stability regions for q vs. r+. The vertically shaded area corre- 
sponds to the regions where AdS4 is the preferred solution, and the horizontally shaded area 
corresponds to the region where the black hole solution is a local minimum of the Helmholtz 
action. 



VI. PHASES OF R-CHARGED BLACK HOLES IN D = 7 

We quote the black hole solution of A^ = 4 gauged supergravity in D = 7 (While the 
solution has not been derived from the actual D = 7 Lagrangian the structure of this solution 
is in close analogy with the structure of the R-charged solution in D = 4 and D = 5.): 

ds' = -{H.H^r^/'fdt' + {H.H^f' l^r'dr' + r'dQ,,k) , (58) 

where 



f = k-^,+g'r'H,H,, iJ, = l + ^, (^ = 1,2). (59) 

2 

We have chosen G7 = ^. We are most interested in the case k = 1. The qi and g^ are again 
related by (6) and the mass is specified as: 

M = f/i + gi + g2 = |« + rlg^.H^) + gi + gs , (60) 

where again in the second equality the radial coordinate at the outer horizon r^ is defined 
as the largest non-negative zero of /(r) (59). 

The parameterization of the thermodynamic quantities in terms of of g^ and r_|_ yields 
the following expressions: 

S = l?r = ^r+Jiri + gi)(r| + g^) , (61) 



AG, 


= 7rr+y'(4 + gi)(r| + g2) , 




^rl^{ri + q,){ri + q,) 



Th 3r^ + 2r^ + 4(gi + g2)-gig2' ^ ' 

The physical charges are then determined in terms of qi and r+ in the following way: 
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gf = qi{ri + qi 



1 



i + -^K + g,; 



, (j 7^ , 



(63) 



and the electric potentials at the outer horizon take the following form: 



Kir^) 



ri + Qi 



1,2 . 



(64) 



In the following we set L 



1. L, which is related to the asymptotic cosmological 



constant A = — -p-, can be restored by replacing r+ by -^ and Qi hj j^. 

A. Grand- Canonical Ensemble and Stability Constraints 

Ig for the black hole solution, when expressed explicitly in terms of qi and r+, takes the 
form: 



tBH 



A 

2rl 



(r+ + gi)(r^ + g2) 



(65) 



where (3 is determined by (62). 

For one non-zero charge only, the thermodynamic stability constraint takes the following 
form: 



3r^ + r% + 2qrX - 2r^ + 3gr^ - g^ > . 
The critical line has two branches: 



(66) 



q±{r+) =rl{rl + \ + ^l^rX + 16r^ + l) , (67) 

and the region of local stability is q-{r^) < q < g+(r_|_). g = has the critical point at 
r+ = w |. Note also that as r_|_ ^ the critical line is pushed to q± -^ 0. 

The phase transition between the black hole solution and the pure AdSf takes place at 
go(^+) = 0, where 



go(r+) = rl-rl . 



(68) 



For q > go(^+) the black hole is the dynamically preferred solution. As r+ ^ the critical 
line is pushed to go ~^ 0. When go = 0, the Hawking-Page transition is at r_|_ = 1. In 
the region g > g+(r+) > go(r_|_), the black hole solutions are unstable, but AdSy is not 
dynamically preferred either. 

In Figure 5 the stability regions are shown in a plot of g vs. r^. 
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FIG. 5. The stability domains as q vs. r+ for the grand-canonical ensemble of D = 7 R-charged 
black holes with one non-zero charge. The vertically shaded areas correspond to the regions where 
the AdSj is the preferred solution and the horizontally shaded area correspond to the regions where 
the black hole solutions are local minimum of the Gibbs action. The checkered area is the domain 
of common overlap. The unshaded area is where neither black holes nor anti-de Sitter space are 
stable. Powers of L can be restored by replacing r-|_ with ^ and q with -X. 

B. Canonical Ensemble and Stability Constraints 

Repeating the same analysis for the canonical ensemble yields the Helmholtz action (23) 
which for the black hole solution takes the following form: 

(3 



tBH 
'■H 



2rl 



— r , + r , 



(3r^ + 4r;)(gi + g2) + 7gig2 



(69) 



with /5 defined in (62). 

In the case on single non-zero charge, the local stability constraint leads to: 
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The critical line has two branches 
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with g_(r+) < q < g+(r+) corresponding to the stable region. Note that as r+ 
critical lines are pushed to g± — > 0. 

The Hawking-Page transition between the black hole solution and the pure AdSj solution 
is now determined by 5'o(^+) = 0, where 
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(72) 



For q > qo{r^) the black hole solution is the dynamically preferred solution. As r+ — > 0, 
qo -^ 0. In Figure 6 these stability domains are shown in a plot of q vs. r_|_. 
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Figure 6 

FIG. 6. The stability domains q vs. vj^ for the canonical ensemble of Z) = 7 R-charged black 
holes with one non-zero charge are given. The vertically shaded areas correspond to the regions 
where the AdSj is the preferred solution and the horizontally shaded area correspond to the regions 
where the black hole solutions are local minimum of the Helmholtz action. The checkered area is 
the domain of common overlap. The units of L = 1 can be restored by replacing r. 
q with ■^. 
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VII. R-CHARGED BLACK HOLES IN D = b {D 

(M5)-BRANES 



7) AND SPINNING M2 



The Kaluza-Klein reduction described in Section IV of spinning D3-branes to R-charged 
black holes in AdS^^ can be repeated for spinning M2- and M5-branes [23], and it leads to 
R-charged black holes in AdS/^ and AdSj. As in Subsection IV B, one can also obtain the 
precise translation from the mass M, charges gj, and the entropy of large R-charged black 
holes in AdS^ [AdSj], to the energy E above extremality, the angular momenta Ji, and the 
entropy S of near-extreme spinning M2-branes [M5-branes] . 

As a consequence of this equivalence, the thermodynamic stability constraints derived 
for the general spinning M2-branes [and M5-branes] should coincide with those of large 
R-charged black holes in D = 4 [and -D = 7], respectively. For the case of only one non- 
zero charge the grand-canonical ensemble (51) [(66)] and canonical ensemble (51) [and (66)] 
stability constraints reduce in the case of large R-charged black holes (r+ ^ 1 and g ^ 1) 
to the following stability constraints: 



3r+ — 2g > , D= 4, grand-canonical 
r+ + 2g > , D= 4, canonical , 



(73) 
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which agree precisely with the respective stability constraints [7,6] for spinning M2-branes 
[and M5-branes] with only one angular momentum turned on. 

Finally, we would like to note a feature of the asymptotics of the stability conditions for 
small and large g which may be only numerological coincidence, but seems to us intriguing. 
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There is an exponent Up which characterizes the scahng of the entropy S of near-extreme 
non-dilatonic p-branes with the number N of branes: S ~ A^"p. For p = 3 (the D3-brane), 
^3 = 2, which is interpreted as evidence that there are indeed on the order of iV^ degrees 
of freedom on the world-volume, as expected in a conformally invariant gauge theory with 
gauge group SU{N). For p = 2 (the M2-brane), 0:2 = 3/2. For p = 5 (the M5-brane), 
0:5 = 3. Cross-sections of minimal scalars falling into any of these branes have the same 
A^"p scaling [24,25]. We note for the record that Up = ^5f where d is one plus the number 
of spatial dimensions perpendicular to the brane world- volume. 

In the grand-canonical ensemble, the ratio of charge to mass where stability is lost for 
large R-charged black holes in D = 4 [Z^ = 7] is most efficiently expressed as .p^ — '^ 



(0:2) ^^"^ l\l ~ y^ ~ ('^5) ^^^]' ^^ ^^^ other hand, in the g ^ limit, the critical mass 

where stability is lost can be determined from r+ = 4= [r+ = w|] in D = A [D = 7]. Thus 
the D = 4 and D = 7 results are in a peculiar numerological sense dual to one another. For 
the case of R-charged black holes in D = 5, the large black hole stability constraint and that 

of g = are self-dual: \/ ^ = "75 ~ ('^s)^^^^ fo^ large g, and r+ = 4^ for small q. 

VIII. CONCLUSIONS 

One of the motivations for this work was the question, "What happens in a quantum 
theory of gravity when a black hole becomes thermodynamically unstable?" String theory 
in principle provides us with a venue to address this question. In general, however, the 
black holes of which one has a good string theory description are thermodynamically stable. 
A possible exception, first noted in [5], is D3-branes with spin, by which we mean angular 
momentum in a plane transverse to the world-volume. These objects, which we have shown 
to reduce to charged black holes in AdS^, can become unstable due to excessively large 
spin — or, in the five- dimensional picture, excessive charge. 

We had initially hoped to show that the charge-driven instability of black holes in AdS^ 
could be understood in terms of the Hawking-Page transition: that is, for the black holes 
which were thermodynamically unstable, there would be a transition to a gas of charged 
particles in pure anti-de Sitter space. To a first approximation, this transition occurs when 
the Euclidean action of the black hole solution rises above the Euclidean action for pure 
anti-de Sitter space subjected to the same boundary conditions at spatial infinity. 

This hope was only partly substantiated, and an intricate phase diagram realized, 
through analytic treatments of the local thermodynamic stability conditions and of the 
criterion for the Hawking-Page transition. Holding temperature and electric potential fixed, 
we found in particular that there is a region of phase space where the black hole solution 
is locally unstable (that is, it is only a saddle point and not a local minimum of the Eu- 
clidean action), but where pure anti-de Sitter space has a still higher action! In this region, 
neither the black hole solution nor the gas of particles in anti-de Sitter space can be the 
preferred state of the system. The only option we can suggest at this point is that the pre- 
ferred configuration is not spherically symmetric with respect to rotations about its center 
of mass in AdS^, despite the fact that it has no angular momentum in AdS^ — spin is angular 
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momentum in the S^ directions. For instance, one might hope for a multi-black hole config- 
uration. Fragmentation of highly charged black holes would be a truly novel phenomenon, 
but it does have some intuition behind it: near-extreme D3-branes which spin too fast seem 
likely to split apart into chunks which move out in the radial direction (which becomes a 
dimension in AdS^ upon Kaluza-Klein reduction). The chunks could carry off some of the 
spin as orbital angular momentum in the S^ directions. The ten-dimensional geometry we 
are envisioning is neither stationary nor static, but it does seem that its five-dimensional 
analog is a multi-black hole geometry. We hope that progress can be made in finding such a 
five-dimensional geometry, or in proving that it cannot exist by some extension of the known 
uniqueness theorems for black hole geometries. 

Before we get carried away with the possibilities for quantum fission of black holes, we 
should point out that the conclusions depend on whether total charge or electric potential 
is held fixed. If total charge is held fixed, then in all the cases we considered {D = 4 and 
D = 7 as well as D = 5), for each black hole geometry which is locally thermodynamically 
unstable, there is a pure anti-de Sitter geometry with the same total charge and periodicity 
of Euclidean time which is entropically preferred in the sense that it has smaller Euclidean 
action. Thus for fixed total charge, our initial hope is realized: black holes which get too 
highly charged to exist in equilibrium with a heat bath undergo a Hawking-Page transition. 
As usual in first order phase transitions, it is also possible for the black hole solution to 
represent a local minimum of the Euclidean action (i.e. a local maximum of the entropy), 
to equilibrate with a heat bath — in short to lead a stable existence — but nevertheless not 
to be the global minimum of the action. If one waits long enough, such "meta-stable" black 
holes should quantum tunnel to a gas of particles in AdS (plus perhaps a locally unstable 
black hole which then decays quickly into the heat bath). 

As we have indicated in Subsection III A, the choice of ensemble depends on physical 
motivation. In the case of spinning branes, where the world-volume is infinite and one is 
really considering a black brane geometry where any part of the world- volume can act as a 
"heat bath" for any other part, the grand canonical ensemble seems the appropriate one for 
judging local stability. For black holes in anti-de Sitter space, the situation is less clear. Our 
gravity intuition is that a conserved charge should be held fixed when studying a localized 
object (the black hole). That is, we should use the canonical ensemble. But if gravity is 
wholly refiected in the gauge theory on the boundary, it seems that, at least in A; -^ O"*" 
limit where the world- volume of the gauge theory is much larger than the cube of its inverse 
temperature, the same arguments that applied to spinning branes tell us that we should 
be using the grand canonical ensemble. We have yet to resolve this issue to our complete 
satisfaction. 

Most of our results were obtained in the case of a single charge non-zero. The number of 
independent charges is the rank of the gauged super gravity's gauge group, or in the language 
of spinning branes the number of mutually orthogonal planes perpendicular to the brane 
world-volume. However, our analysis lays the groundwork for an exploration of multiple 
non-zero angular momenta. In particular, we have obtained explicit closed forms for the 
Euclidean actions, both with charges held fixed and with potentials held fixed, for arbitrary 
numbers of charges up to the maximum number allowed: three for the D3-brane, four for 
the M2-brane, and two for the M5-brane. Using these actions to locate the Hawking-Page 
transition is straightforward: the way we have set up the calculations, the transition takes 
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place when the action is zero. Evaluating the local stability constraints may be feasible in 
the general case for the canonical ensemble, but it is exceedingly calculationally burdensome 
in the grand-canonical ensemble because determinants of matrices with up to five rows and 
columns are involved. In the case of large black holes we have been able to evaluate these 
determinants. This and other issues will be reported on in [6]. 



Note added. As this paper was being completed we received the preprint [26], which 

overlaps with some of our results. 
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